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Abstract 

Abstract We study the evolution of the energy (mode-power) distribution for 
a class of randomly perturbed Hamiltonian partial differential equations and derive 
master equations for the dynamics of the expected power in the discrete modes. In 
the case where the unperturbed dynamics has only discrete frequencies (finitely or 
infinitely many) the mode-power distribution is governed by an equation of discrete 
diffusion type for times of order 0(e~ 2 ). Here e denotes the size of the random 
perturbation. If the unperturbed system has discrete and continuous spectrum the 
mode-power distribution is governed by an equation of discrete diffusion-damping 
type for times of order 0(e~ 2 ). The methods involve an extension of the authors' 
work on deterministic periodic and almost periodic perturbations, and yield new 
results which complement results of others, derived by probabilistic methods. 
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1 Introduction 

The evolution of an arbitrary initial condition of linear autonomous Hamiltonian partial 
differential equation (Schrodinger equation), 



where H is self-adjoint operator, can be studied by decomposing the initial state in terms of 
the eigenstates (bound modes) and generalized eigenstates (radiation or continuum modes) 
of H . The mode amplitudes evolve independently according to a system of decoupled 
ordinary differential equations and the energy or power in each mode, the square of the 
mode amplitude, is independent of time. If the system (JTJ) is perturbed 



where W(t) respects the Hamiltonian structure (W* = W), then the system of ordinary 
differential equations typically becomes an infinite coupled system of equations, so-called 
coupled mode equations. If W(t) has general time-dependence (periodic, almost periodic, 
random,...), the solutions of the coupled mode equations can exhibit very complex behavior. 
Of fundamental importance is the question how the mode-powers evolve with t. Kinetic 
equations, which govern their evolution are called master equations [?], [?] and go back to 
the work of Pauli [?]. A general approach to stochastic systems is presented in [?, ?, ?, ?]; 
see also [?, ?, ?]. Master equations have been derived in many contexts in statistical 
mechanics, ocean acoustics and optical wave-propagation in waveguides. 

We present a theory of power evolution for (J2J), for a class of perturbations, W(t), which 
are random in t. Our theory handles the case where H has spectrum consisting of bound 
states (finitely or infinitely many discrete eigenvalues) and radiation modes (continuous 
spectrum). It is a natural extension of the analysis in our work on deterministic periodic, 
almost periodic and nonlinear systems; see, for example, [?, ?, ?, ?]. Our approach is 
complementary to the probablistic approach of [?, ?, ?, ?, ?]. The model we consider is 
well-suited to the study of the effects of an "engineered" perturbation of the system, e.g. a 
prescribed train of light pulses incident on an atomic system, or prescribed distribution of 
defects encountered by waves propagating along a waveguide; see below. We also give very 
detailed information on the energy transfer between the subsystems governed by discrete 
"oscillators" and continuum "radiation field" . 



id t (j) = H <f>, 



(1) 



idt<j> = (H + W(t))<f>, 



(2) 
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In particular, we study the problem 

idt<p = ( H + eg(t)P ) 0, 



(3) 



where e is small, and H and (3 are self-adjoint operator on the Hilbert space T~C. Hq is 
assumed to support finitely or infinitely many bound states. For example, Hq = — A + V(x), 
where V is smooth and sufficiently rapidly decaying as \x\ — > oo. /3 is assumed to be 
bounded. g(t) is a real valued function of the form of a sequence of short-lived perturbations 
or "defects"; see figure 1. Our methods can treat the case of more general perturbations, 
e.g. W(t,x) = /3(t,x), but to simplify the presentation we consider the separable case 
W(t) = g(t)P(x). 

Models of the above type arise natural in many contexts. Among them are the inter- 
action between an atom and a train of light pulses [?, and references therein] , a field of 
great current interest in the control of quantum systems. Such trains of localized pertur- 
bations also model sequences of localized defects along waveguides, see [?], [?], introduced 
by accident or design. 

We construct g(t) as follows. Start with go(t), a fixed real- valued function with support 
contained in the interval [0, T] and let {dj}j>o be a nonnegative sequence. Define 



git) = J2 9o{t 



t n ), where 



n=0 



dn 



-- (d + T) + (d 1 +T) 
denotes the onset of the n th defect. 



. . . + (d„_i + T) + d n ,n > 1 



(4) 



(5) 



g(t) 



do 



ga(t-d ) g (t ~ T - d - di) 90 (t - (n - 1)T - d ■ 



rfl 






Figure 1: Train of short lived perturbations or "defects". The onset time for the n defect, 
t nj is given by ©. 

Note that, if the sequence {4/}j>o is periodic then g(t) is periodic. In this case, the 
system Q has already been analyzed by time-independent methods [?] or, more recently 
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and under less restrictive hypothesis, in [?, ?]. For {dj}j> quasiperiodic or almost periodic 
(see [?, ?] for a definition) the situation is more delicate. In [?] we treat a general class of 
almost periodic perturbations of the form: 

W(t) = ^cos(/^)&, (6) 

i>o 

with appropriate "small denominator" hypotheses on the frequencies {/ij}. We leave it 
for a future paper [?] to consider the case of almost periodic {dj}j> and to explore the 
connection with the results in [?]. We note that a particular case has already been treated 
in [?, Appendix E]. 

Note that in [?] and [?] the numbers do, d\, . . . , are equal to a fixed constant and go(t) 
is random while in our model do,di,..., are random and go(t) is fixed. This is another 
sense, in which our results complement those in the existing literature. 

The paper is divided in two parts. The first part treats stochastic perturbations of 
Hamiltonian systems with discrete frequencies and then second part extends these results 
to the case where the unperturbed system has discrete and continuous frequencies. The 
stochastic perturbation is of order e and then the vector P(r) G i 1 , whose components 
are the expected values of the squared discrete mode amplitudes (mode-powers), satisfies 
on time scales t = 0(e~ 2 ) or equivalently r = 0(1), the master equations of diffusion or 
diffusion-damping type. Specifically, if H has only discrete spectrum (finite or infinite) 
then 

d T P(r) = -BP(t), B>0 (7) 
which has the character of a discrete diffusion equation, i.e. 

E P *( r ) = E P *(°)' fP'P = -(P,BP) <0. (8) 
k k T 

If H has both discrete and continuous spectra, then 

d T P(r) = (-B - T)P(r), B > 0, T = diag( 7fe ) > (9) 

for which 

J> fe (r) < e"^^P fc (0), (10) 

k k 

where 7 = min^ 7^. 

In sections 2 and 3 we study (jSJ) under the hypothesis that H has no continuous 
spectrum (i.e. no radiation modes) and in section 4 we generalize to the case where H 
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has discrete and continuous spectrum. In section |21 we present the main hypotheses on H 
and g (t) and study the effect of a single short lived perturbation. In section El we present 
our hypotheses on d ,di, . . . , and analyze the effect of a train of perturbations (J3JIIJ). We 
show that if do,di, . . . , are independent random variables with certain distributions, see 
Hypothesis (H4) and Examples 1 and 2, diffusion occurs in the expected value for the 
powers of the modes. Specifically, if we start with energy in one mode, then, on a time 
scale of order 1/e 2 , one can expect the energy to be distributed among all the modes. 
In section 0] we analyze equation (jSJ) under the hypothesis that H has both discrete and 
continuous spectrum (i.e. supports both bound modes and radiation modes). We prove a 
result similar to the nonradiative case but now bound state- wave resonances lead to loss 
of power. The effect of our randomly distributed deterministic perturbation is very similar 
to the one induced by purely stochastic perturbations, see [?, ?, ?], but quite different from 
the effects of time almost periodic perturbations, see [?, ?]. Section |3] is dedicated to such 
comparisons. 
Notation 

1) (x) = y/TTx 1 

2) Fourier Transform: 



3) We write ( + c.c. to mean ( + (, where ( denotes the complex conjugate of (. 

4) w' denotes the transpose of w. 

5) [q\ denotes the integer part of q. 

Acknowledgement: We would like to thank G. C. Papanicolaou and J. L. Lebowitz for 
helpful discussions concerning this work. E.K. was supported in part by the ASCI Flash 
Center at the University of Chicago. 

2 Short lived perturbation of a system with discrete 
frequencies 

In this section we consider the perturbed dynamical system 



where H has only discrete spectrum and go(t) is a short-lived (compactly supported) 
function. We study the effect of this perturbation on the distribution of energy among the 




(11) 



idt<t>(t) = H Q <t>(t) + eg {t)/3<i>(t,x), 



(12) 
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modes of H . Here and in section 0] we are extending the results in [?] to multiple bound 
states but under an additional "incoherence" assumption; see (JTHJ). 
Hypotheses on H , (3 and go(t) 

(HI) H is a self adjoint operator on a Hilbert space TC. It has a pure point spectrum 
formed by the eigenvalues : {Aj}j>i with a complete set of orthonormal eigenvectors: 
{ipj}j>i ■ 

Hoipj = Xjtpj , = 5ij (13) 

(H2) (3 is a bounded self adjoint operator on Ti and satisfies ||/3|| = 1. 

(H3) g (t) G L 2 (R) is real valued, has compact support contained in [0,T] on the 
positive real line and its L 1 -norm, denoted by ||<7o||i is 1- Thus its Fourier transform has 
L°°-norm bounded by 1. 

Note that one can always take \\j3\\ = 1 and \\go\\i = 1 by setting e = \\go\\i • thus 
incorporating the size of go/3 in e. Therefore, under assumptions (H2-H3), e in (|12|) measures 
the actual size of the perturbation in the L l (R,7i) norm. Our results are perturbative in 
e and are valid for e sufficiently small. 

By the standard contraction method one can show that (|T2*j) has an unique solution 
4>{t) G H for all t G 1. Moreover, because both H and go(t)/3 are self adjoint operators, 
we have for all t G R : 

urn = wmi (14) 

We can write (j>{t) as a sum of projections onto the complete set of orthonormal eigenvectors 
of H : 

0(t, x) = a i( t )V'i( a; )' ( 15 ) 

3 

By Parseval's relation 

EM*)I 8 = II^)II 2 = II^°)II 2 ( 16 ) 

3 

Now (|12jl can be rewritten as 

id t a k {t) = X k a k {t)+ego{t)Y, a 3( t )(^k,^ j ), fee {1,2...} (17) 

3 

where (-, ■) denotes the scalar product in TC. 

Hence the equation fT2|) is equivalent to a weakly coupled linear system in the ampli- 
tudes: ai, a 2 , ■ ■ ■ , (JUj). 
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Since the perturbation size is e we expect, in general, that the change in energy in 
the k th mode, |afc(t)| 2 — |afc(0)| 2 , to be of order e. However with a suitable random initial 
condition we can prove more subtle behavior. 

Suppose that there exists an averaging procedure applicable to the amplitudes: a±, a 2 , . . . 
of the solutions of (fT2|) . denoted by 

a(t) i — ► E(a(t)) G C. 

We now state a fundamental result, applied throughout this paper, for a single defect 
which is compactly supported in time: 

Theorem 2.1. Assume the conditions (H1)-(H3) hold and the initial values for Alty) are 

such that 

E (aj(0)a fc (0)) = whenever j k. (18) 
Then for all t > sup{s G R | go(s) ^ 0} and k G {1,2,.. .} we have 

P k (t) - P k (0) = z 2 J2 l^il 2 |^o(-A fei )| 2 (P,(0) - P k (0)) + 0(e 3 ), (19) 

3 

where 

P k (t)=E(\a k (t)\ 2 ) 

denotes the average power in the k th -mode at time t, a k j = (ip k , ftipj), go denotes the Fourier 
transform of go and A k j = X k — Xj. 

Note that (|T9"j) can be written in the form: 

P k (t) =T £ P k (0) + O(e 3 ), (20) 

where 

T £ = I-e 2 B; B>0 (21) 
I is the identity operator (matrix) and B is given by 

U - (h \ h _ / -Kil 2 |£o(-A fcj )| 2 , for j^k, 



J2i,i^k \ a ki\ \9o(-^h)\ , for j = k 
In section El we will discuss and use the properties of B and T £ . 

Proof of Theorem 12.11 In the amplitude system, (|17p. we remove the fast oscillations by 
letting 

a k {t) = e~ iXkt A k {t), (23) 
7 



Note that by (fTKjl 



Now (fTTjl becomes 



where 



£|^(i)| 2 = ||«K0)|| 2 (24) 
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id t A k (t) = eg (t) akje^Mt), (25) 



Akj = Afc — Xj, (26) 
ctfcj = (ipk, = oijk- (27) 

The above system leads to the following one in product of amplitudes, A k (t)Ai(t): 

d t (A k (t)Mt)) = ieg (t)J2 a Ji eiAj,tA k(t)Mt) 

j 

iegoit^ockje^A^Mt), (28) 
i 

In the particular case k = I we have the power equation for each mode: 

d t \A k (t)\ 2 = ieg (t) ajke^A^Ajit) + c.c. . (29) 
j 

Note that the sum in (|29|) commutes with time integral and expected value operators. This 
is due to (j2*4^) and the dominant convergence theorem, see for example [?] . Indeed consider 

m 

f m (t) = J2 a jke^ kt A k (t)A J (t)g (t). 

3=1 

From (jl5j) we have for all t G R 

lim f m (t) = ((f)(t), f3ip k )ak(t)g (t). 

From (J2~^ and the Cauchy-Schwarz inequality \(a, b)\ < ||a|| ||6||, we have for all t G R 

\f m (t)\ < H(0)\\ 2 \9o(t)\. (30) 

The right hand side of (J3*U|) is integrable and the dominant convergence theorem applies. 
A similar argument is valid for expected values. Therefore, from now on, we are going to 
commute both time integrals and expected values with summations like the one in (}29|) . 
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We integrate from to t > sup{s G K. | <?o( s ) ^ 0} and integrate by parts the right 
hand side. The result is: 

\A k (t)\ 2 -\A k ($)\ 2 = ieJ2a jk [ ' g (s)e l ^ s A k (s)A,(s) + c.c. 

Jo 

/do 
g^e^drA^A^s) |^ + c.c. (31) 
j 

pt pOO 

+ kVa jt / g (T)e tA ^ T dTd s (A k A j )(s)ds + c.c. . 

The boundary terms are 

^ (r)e iA ^ r drA fe (s)Z i (s) |^ + c.c. 

= ie^a ife y o (-A ijfc )A fc (0)Z i (0) + c.c., (32) 
j 

where go denotes the Fourier Transform of go; see (fTTj) . Note that upon taking the average, 
using (fTHj) and the fact that <?o(0) is real, these boundary terms vanish. 
Into the last term in (pHj) we substitute (J2~%j) : 

pt poo 

ieJ2 a jk / g (r)e iA ^ T dTd s (A k A j )(s)ds = 

pt poo 

= +\e\ 2 y^a jk a kp j I g (T)e iA ^ T dTg (s)e iA ^ s A p (s)A j (s)ds 

pt poo 

- | £ | 2 Va Aj / / g (r)e^ T drg (s)e^ s A k (s)A q (s)ds. (33) 
We again integrate by parts both terms in ()33|) : 

pt poo 

ieS"a 3k / g (T)e iA ^ T dTd s (A k A J )(s)ds = 

j JO Js 

poo poo 

- \e\ 2 ^a jk a kp g (s)e iAk ? s g (r)e lA ^ 7 drdsA^ujA^u)^ 



1,P 

roc 



I u=t 
^qK^J \u=0 



poo poo 

+ \e\ 2 ^a jk a q] / go(s)e iA ^ s / g (r)e iA ^ T drdsA k (u)A q 

moo 
g (T)e iA ^ T dTg (s)e lAk ^dsd u (A p A 3 ) [u)du 
hp 

moo 
go(T)e tA ^ T drg (s)e iA ^ s dsd u (A k A q ) (u)du. (34) 



Note that the boundary terms calculated at "it = t" are zero since t > sup{s £ 
^ I 9o( s ) 7^ 0}. Upon taking the expected value and using (JTSj) the only boundary terms 
contributing are the ones for which u = and j = p in the second row of 



/*00 POO 

/ / 9o(r)e iA ^ T drg (s)e iA ^dsE (^(O)! 2 ) + c.c. 

3 J ° Js 

poo poo 

= V|a fe /E(|A,-(0)| 2 )-23J / / g (T)e iA ^drg (s)e iA ^ds (35) 

• JO Js 

and the ones for which u = and q = k in the third row of (|34|): 

/*oo /»oo 

S Kf / / go(r)e lA ^dTg (s)e lA ^ s dsE {\A k (0)\ 2 ) + c.c. 

3 J ° Js 

/OO POO 
J go(T)e lA ^ T drg (s)c lA ^ds (36) 



To compute (j35H36|) we use the lemma: 

Lemma 2.1. If go(t),t £ R is square integrable with compact support included in the 
positive real line then for all A £ R the following identity holds 

/oo poo 
J g (r)e aT drg (s)e- iXs ds = \g (-X)\ 2 . 

Proof. For any A £ R we have: 

g (r)e lXT drg (s)e- iXs ds 



Js 



poo poo 

= Hm / g (T)e l ^ T dTg (s)e~ as ds 

£ \°Jo Js 
-i poo poo poo 

= — hm / g (s)e- lXs ds / g (») / e^+^drdfi 
= — hm [°° g (s)e- lXs ds f°° ^ e^ +A+ ^/i 

= -Llixn f°MH-ti d 
2ne\oJ_ 00 fi + X + te 

= i| 9 o(-A)| 2 + fp.V.fM^f^ (37) 

The last relation in (J37)l is the Plemelj-Sohotsky's formula for (temperate) distributions: 

lim — = P.V. iir5(x) 



e\o x + ie x x + iO 
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Note that <7o(/f) £ C°°(R) H L 2 (R). Since ()37|) is already decomposed in its real and imagi- 
nary part the lemma follows. [] 

Into the triple integral terms of (}3"4*j) we again substitute ([28)1 . Then one can show that 
the 1-norm of this correction vector is dominated by \e\ 3 ||<7o||i ||/5|| 3 ||0(O)|| 2 . Hence, it is 
of order C(|£:| 3 ). 

Thus, after applying Lemma l2~Tl to (j35H36|) and using (|3*T|) we arrive at the conclusion 
of Theorem 12.11 j 

3 Diffusion of power in discrete frequency (nonradia- 
tive) systems 

In the previous section we calculated the effect of a single defect on the the mode-power 
distribution. In this section we show how to apply this result to prove diffusion of power 
for the perturbed Hamiltonian system, (j2J), where g(t) is a random function of the form (j3J), 
defined in terms of a random sequence {dj}j>o. In particular, the sequence {<i,}j>o will be 
taken to be generated by independent, identically distributed random variables. This will 
be result in a mixing the phases of the complex mode amplitudes, after each defect. 

We assume that (H1-H3) are satisfied. The following hypothesis ensures that (JTBJ) 
holds before each defect, thus enabling repeated application of Theorem 12. II 

(H4) d , di, ... are independent identically distributed random variables taking only 
nonnegative values and such that for any I G {0, 1, . . .} and j ^ k G {1, 2 . . .} we have 

E ( e i(A ^~ Afe)di ) = 

where E(-) denotes the expected value. 

Clearly (H4) requires the eigenvalues to be distinct but aside from these we claim that 
for any finitely many, distinct eigenvalues Ai,A2, ...,A m there exist a random variable 
satisfying (H4). 

Example 1 (finitely many bound states) Given Ai, A2, . . . , A m distinct choose the 
random variables di, I — 0, 1, . . . to be identically distributed with distribution d : 

d= d i k 

l<j<k<m 

where djk are independent random variables such that the distribution of djk is uniform on 
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the interval [0, 2n/\\j — \ k \]. In this case, for any f ^ k' G {1, 2, . . .} 



e (y^'~ Afc '^) = e ( e *( A j'- A fe')'ijfc j 

\l<J<fc / 



a<i< 

= Y\_ E(e i( V~ A fc'Kfc) 

i<i<fc 
= 

since E (e^'-VH'fc') = o. 

Another choice is to consider discrete cL^'s. Namely, take djk to be the discrete random 
variable taking each of the values and 7r/|Aj — A&I with probability 1/2. A concrete 
example is, in the case we have three eigenvalues Ai < A2 < A3, to choose d to be the 
random variable taking each of the eight values: 

7T 71 7T 



A2 — Ai A3 — Ai A3 — A2 



A2 — Ai A3 — Ai J VA2 — Ai A3 — A2 / — ^1 ^3 — ^ 

1 1 1 

71 I " — + " — + 



A2 — Ai A3 — Ai A3 — A2 / 
with probability 1/8. 

(H4) does not restrict us to system with finitely many bound states: 
Example 2 (infinitely many bound states) Let the quantum harmonic oscillator in 
one dimension: 

H = -H-d 2 x +u; 2 x 2 , 16R, 
2 

be the unperturbed Hamiltonian. Then 

A n = huj{n+ 1/2), n = 0,1,2,..., 

see for example [?]. Note that (H4) holds provided that we choose di, 1 = 0,1,... to be 
identically and uniformly distributed on the interval [0, 27i/(huj)]. 

Note on degenerate eigenvalues: As discussed above (H4) cannot be satisfied in the case 
Ho admits degenerate eigenvalues. However, at least in some cases, our theory can be 
applied. In general the degeneracy is a consequence of the symmetries of Hq, i.e. the 
existence of a self-adjoint operator, say L, commuting with Hq, [L, Hq] = 0. To recover our 
results it is sufficient to assume that (3, the "space-like" part of the perturbation, respects 
the symmetry, i.e. commutes with L. One can now factor out L, i.e. work on the invariant 
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subspaces of L where H is nondegenerate. Along the lines of Example 2 one can consider 
the quantum harmonic oscillator in three dimensions which has a spherically symmetric 
Hamiltonian and degenerate eigenvalues, see for example [?] . If /? is spherically symmetric 
then it only couples bound states with the same angular momentum. Hence the problem 
reduces to subsystems consisting of bound states with the same angular momentum but 
different energy, therefore nondegenerate. The choice we made in Example 2 will satisfy 
(H4) in each of the subsystems. 

3.1 Power diffusion after a fixed (large) number of defects 

Theorem 3.1. Consider equation h!2\) with g of the form O). Assume (H1-H4) hold. 

Then the expected value of the power vector after passing a fixed number of perturbations 
"n" satisfies 

p(n) =T e n P(O)+0(n£ 3 ), (38) 

where T £ is given in 1121]) 

P fc (n) =E(Mt)| 2 ), A; = 1,2,... (39) 
t„,_i + T < t < t n , (t ranging between the n th and (n + l) st defects 

Proof. We will prove the theorem by induction on n > 0, the number of defects traversed. 
For n = the assertion is obvious. Suppose now that for n > we have 

P (n) = T «p(o) + 0{ne 3 ). (40) 

We will show 

P (n+i) = T n+ip( ) + O ((n + l)e 3 ) (41) 

by applying Theorem 12.11 to (}4*0|) . In order to apply Theorem 12.11 we need to verify that 
(|18j) is satisfied before the n + 1 st defect. Specifically, we must verify that for any pair k ^ j 

E (a k (t n+1 )a 3 (t n+1 )) = E (a* (nT + £3K d *) % (nT + = 0. (42) 

Using the fact that d n+ i is independent of do + d\ + . . . + d n , and (H4) we have: 

E (a k aj (nT + = E ( a ^j (nT + ££ =0 4) e*^-^)^ 1 ) 

= E (auaj (nT + ELo^)) E (e^^^ 1 ) = 0. 
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Thus (J42j) holds and all the hypothesis of Theorem 12.11 are now satisfied. By applying it 
and using (j4~Uj) we have 

p (n+i) = T e P {n) + 0(e*) 

= T e (T"P(0) + 0{ne 3 )) + 0{e 3 ) 
= T? +1 P(0) + O((n + l)s 3 ). 

Hence (f4*Uj) implies (}4~Tj) . This concludes the induction step and the proof of Theorem 13.11 
is now complete. 

In the next two Corollaries we describe the asymptotic behavior of the vector of expected 
powers when the number of defects n tends to infinity. Note that after a possible reordering 
of the eigenvectors ipi, tpz, . . . , of H , the operator B given by (J2~2"j) might look like 1 : 

B = di a g[B 1 ,B 2 ,...,B g ,...}, (43) 

where B\, E>2, . . . , B q , . . . are square matrices (linear operators) of dimensions rni, 777,2, . . . , m q , . . . , 
1 < m q < oo, q = 1, 2, .... In linear algebra terms this means that B is reducible. In terms 
of the dynamical system (|3*Hj) generated by T £ = I — e 2 B it means that, after a possible 
reordering, the first vii\ bound states of Hq are isolated from the rest. The same is valid 
for the next rri2 bound states, etc. To understand the evolution of the full system it is 
sufficient to analyze each of the isolated subsystems separately. They all evolve according 
to (J3~%j) with T e = I — e 2 B q and B q given by (f2~2"|) but the indices span only a subset of the 
eigenvectors ipi, ip2, ■ ■ ■ of Hq. The main difference is that now B q is irreducible. In what 
follows we are focusing on one such subsystem and drop the index q. 

Corollary 3.1. If the subsystem has a finite number of bound states, say m, then 

( P(0), ifn<^e~ 2 
lim P {n) = { e- Br P(0) if n = re~ 2 , (44) 

I 1,-- .,1)', if e- 2 «n«\e\-* 

where E = -Pi(O) + ^(0) + . . . + -P m (0) is the expected total power in the subsystem and it 
is conserved. 

Proof. We use the following properties of the irreducible matrix B: 

(Bl) B is self adjoint and B > 0; 

1 For such a decomposition to occur it is sufficient that Hq and (3 have common invariant subspaces 
Tii C 7i, 7^2 C 7i, . . . , Ti.q C 7i, . . . 
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(B2) is a simple eigenvalue for B with corresponding normalized eigenvector 

r = -U(l,l,...,l)'. (45) 

These properties are proved in the Appendix. 

Let (3q = 0, Pi, (3%, . . . , P m -i be the eigenvalues of B counting multiplicity, and let 
r , 7"i, . . . , r m _i be the corresponding orthonormalized eigenvectors. By (Bl) and (B2) 
Pi, Pi, ■ ■ ■ , Pm-i are strictly positive. Let 

R = [7*0,7"!, . . . ,r m _i] 

be the matrix whose columns are orthonormalized eigenvectors of B and let R! be its 
transpose. Then 

R'BR = diag[Ab A, An-i] 
R'R = I = RR'. 



It follows that 

T n = (l-e 2 B) n = R[R'(l-E 2 B)R] n R' 

= tfdiag [(1 - e 2 p T, (1 - e 2 Pi) n , . . . , (1 - e 2 p m _i) n ] R'. 

We now study lim^oo T™ for the three asymptotic regimes of (pHJ) . Note that for < k < 
m — 1 we have: 



Consequently, 



lim (l~£ 2 Pk) n 



lim (1 - e 2 P 



2n \n 



n— roo,e z -n=r 



lim (1 - e z p 



2n \n 



n— >oo,e^n— >oo 



lim {l-e 2 p k ) n 



n— »oo,£^n— >oo 



i2diag[l, !,...,!]# = 1 



lim T £ n = I #diag [e~^ T ,e 



-An-1T 



0, Pk>0 

1, p k = o 



,-Bt 



if n < £ 2 
if n = tc~ 2 



i?diag[l, 0, 0, ... , 0]R' = projection onto r if e <C n <C |e| 

where r is defined in (j43|) . 

Substitution of (|4T|) into (JHHj) completes the proof of Corollary 13.11 | j 



(46) 



(47) 
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Corollary 3.2. // the subsystem has an infinite number of bound states, then 

limPW = ( P i B °)' .f n<<£ :l , (48) 
rwo [ e ^ r P(0) if n = re z y ' 

For n ^> e~ 2 the limit in £ 2 is 0, while the limit in £ l does not exist. More precisely, 
although the total power in the subsystem is conserved, 

oo 

Y,Pi n) = E, Wn>0, (49) 
fc=i 

{P (n) } does not converge in t due to an energy transfer to the high modes. In particular, 
for any fixed N > 1 : 

oo 

lim V = E, 

n^oo k 
k=N 

N 

limVP fc (n) = 0. (50) 

n— »oo ' * 
k=l 

We note that similar results have been obtained in [?] but for different types of random 
perturbation. 

Corollaries 13.11 and 13.21 show that, on time scales of order 1/e 2 , the dynamical system 
is equivalent with 

8 t P(t) = -BP(t). (51) 

Moreover the definition of — B in (J22|) together with — B < and e~ B unitary on i 1 implies 
that the flow (|51|) is very much like that of a discrete heat or diffusion equation. 

In conclusion the number of defects encountered should be comparable with 1/e 2 to 
have a significant effect. Once they are numerous enough, the defects diffuse the power in 
the system. If the number of defects is much larger than 1/e 2 the power becomes uniformly 
distributed among the bound states. 

Remark 3.1. Hyptothesis (H4) is important. If we do not assume (H4) then the cor- 
rection term for each defect is of size e, since the boundary terms \3ty no longer vanish. 
Consequently the correction term in the main result \3ty is 0(ne) which on the "diffusion 
time scale" n ~ e~ 2 is very large. 

Proof of Corollary 13. 21 In the case of an infinite number of bound states B has the following 
properties, see the Appendix: 
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(Bloc) B is a nonnegative, bounded self adjoint operator on £ 2 with spectral radius less or 
equal to 2; 

(B2oo) is not an eigenvalue for B; 

(B3oo) B is a bounded operator on I 1 with norm ||-B||i < 2; 

(B4JX,) For |e| < 1 the operator T e = (I — e 2 B) transforms positive vectors (i.e. all compo- 
nents positive) into positive vectors and conserves their I 1 norm. 

We are going to focus first on £ 2 results. Based on the spectral representation theorem, 
see [?], we have for any Borel measurable real function / : 

f(B)= f f(s)dn(s). (52) 
Jo 

Here d^(s) is the spectral measure induced by B. Note that B2oo implies the continuity of 
fx(s) at zero. 
Now 

T: = (I — e 2 B) n = f (1 - e 2 s) n d^s) 
Jo 

and 2 

limT £ n = lim/ (l-e 2 s) n dfx{s)= [ lim (l - e 2 s) n d/i{s). (53) 

n-»oo e-O/d J q e^O 

For the last equality we used the dominant convergence theorem with |1 — e 2 s\ n < 1 for 
< s < 2, |e| < 1 and J Q 2 ld/i(s) = I. Using (jUJ), with s replacing (3k, we have that (}53j) 
becomes 

( J 2 ldfi{s) = 1 if n < e 2 

limT e n =^ f* e - aT dn(s) = e~ Br if n = re" 2 , (54) 

[ /x(0+) - ^(0) = if e~ 2 < n < |e|~ 3 

where we used (J52*Jl and the continuity of fi(s) at zero . 
Plugging (jHljl in (jHEJ) gives the required results in £ 2 . 
For the results in £ we use series expansions: 



(I - e 2 B) n = I + ( i ) ^(-5) + ( 2 ) ^(- 5 ) 2 + •••+( n ) 



(55) 



Since < 2, (see property B3oo), the finite series above is dominated in £ x operator 

norm by: 

2 / 71 \ i /o_2\2 ( n \ , , /r>_2\n I n \ /i , r. 2\n / „2ne 



1 + 2e^ x j + (2^ ^ 'J J+... + (2s 2 ) n ^ J = (1 + 2e 2 ) n < e 2n£ . (56) 
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As n — > oo the series in (J56j) becomes infinite. However, as long as n < r/e 2 , r > fixed, 
the sum in f!56|) is finite and hence that in ()55|) is convergent. Now for each k = 1,2,... 
the (k + l) st term in the series f!55|) has the property: 

( n \ 2k( , k _ j if 
J™ U J ( } ~ I i(S) k if n = re- 2 

Hence by the Weierstrass criterion for absolutely convergent series we have: 

,. r (J 2Rr / 1-0 + 0-... = ! ifrKe- 1 

lim T' — hm I — e B) = < (l . m 3 D „ „ 57 

v ' ' I- tB + - + . . . = e~ rB iin = TE- 2 



n^oo 



It remains to prove that as n — » oo,e 2 n — > cxd, {P^} does not converge in £ l . Let 
p(o) ^ £i p| £2 d eno ^ e a vector with positive components, and consider the sequence: 

p(n) = T np(0) g p, £ 2_ ( 58 ) 

By the third part of (|54j) . ||P^ n ^||2 — ^ 0. Assume now that there exists P <E i 1 such that 
Upto-Pld = 0. Since both i 1 and£ 2 convergence imply convergence of each component, we 
deduce that P = 0. On the other hand, by P^ = T e P^ n_1 \ n — 1, 2, . . . and property 
B4oo, we deduce that P (ri) is a positive vector for which ||P (n) ||i = ||-P (0) ||i = E > for all 
n > 0. Consequently P is a nonnegative vector with ||P||i = E > 0, a contradiction. The 
proof of the Corollary is now complete. [] . 

3.2 Power diffusion after a fixed (large) time interval and a ran- 
dom number of defects 

As pointed out in its statement, Theorem 13. II is valid when one measures the power vector 
after a fixed number of defects "n" regardless of the realizations of the random variables. 
That is after each realization of do,dx,... the power vector is measured in between the 
n th and the (n + l) st defect. Averaging the measurements over all the realizations of 
UQ) Oi, U2j ■ ■ ■ gives the result of Theorem 13.11 What happens if one chooses to measure the 
power vector at a fixed time "t" (i.e. a fixed distance along the fiber)? The answer is given 
by the next theorem: 

Theorem 3.2. Consider equation hlty with g of the form Q). Assume that (H1-H4) 
are satisfied and that all random variables do, di, . . . , have finite mean, variance and third 
momentum. Fix a time t, < t 1 / 1 e: | 3 . Then the expected value of the power vector at 
a fixed time P(t) satisfies 

P(t) = T e n P(0) + C(max{te 3 , e 4/5 }), (59) 
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where n = \t/(T + M)\ denotes the integer part oft/(T + M), T is the common time span 
of the defects and M is the mean of the identically distributed random variables do, di, . . . . 

Corollary 3.3. In this setting, the conclusions of Corollaries \S. 1\ HOI hold with n replaced 
by t. 

Proof of Theorem 13.21 As before, let pV°' be the expected power vector after exactly 
defects. Denote by N the random variable counting the number of "defects" up until time 
t, i.e. 

(N - 1)T + d + . . . + d N ^ <t<NT + d + ... + d N . (60) 
and let 5(e) denote the integer, which grows as e decreases: 

5 = max { 1.39 ( 0/ J — ] e" 6/5 , — ^— Jn log (g' 2 ) + ( -=^-^\ log (e~ 2 ) 1 
1 \a 2 [T + M)J 'T + M v 6V ; \T + MJ &v J j 

S = [5\+l, (61) 

where M, a 2 , respectively p are the mean, variance and the centered third momentum, of 
the identically distributed variables ofo, di, c?2, • • • , and n is the integer part of t/(T + M). 
Note that for t ~ e~ 3 or smaller 5 e~ 2 . The choice of 5(e) is explained below. 
The proof consists of three stages: 

1. P(t) = P( n+ V + 0(e) + 0(5e 2 ) 

2. P("+ 5 ) = + 0(5e 2 ) 

3. PW = T?P(0) + 0(ne 3 ) 

where n = \ t/(T + M)J. The last stage is simply Theorem 13.11 

For the second stage one applies again the previous theorem to get: 

P (n+8) = t^P^ + Q(5e 3 ). 

Now T £ = I — 0(e 2 ) and since 5 e~ 2 stage two follows. 

The first stage is the trickiest. Without loss of generality we can assume that t/(T + M) 
is an integer. Indeed, for n = [t/(T + M)\ we have 

P(t) - P(n(T + d)) = (e(T + M)) = 0(e), 

an error which is already accounted for in this stage. 
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Suppose first n — 5 < N < n + 5, i.e. we condition the expected values to the realization 
of \N — S\ < 0. Then the difference between the conditional expected values of the power 
vector at time t and after n + 5 defects is of order 0(e) + 0(5e 2 ). This follows from the 
fact that the condition n — 5 < N < n + 5 restricts only the realizations of do, d±, . . . djv 
leaving the realizations of dw+i, ■ ■ ■ d n+ s arbitrary; see (J60|) . Hence, as in stage two, the 
conditional expected values satisfy: 

p{n+5) = p{N+l) +(D ^ e 2y 

In addition 

P (N+i) = p(t) + o(e), 

since there are at most 2 defects of size e from "t" up until after the (N + l) th defect. 
Let p{t) denote the power vector 

P {t) = (\a 1 {t)\\\a 2 {t)\\..). 

Recall that by definition P(t) = K(p(t)) and the total power in the system (|T2|) is conserved, 
i.e. 

Ib(t)||i = 5]|a^)| 2 ^|b(0)||i, teR. 

k 

Moreover, 

P(t) = E(p(t) : \N -n\ <5)+E(p(t) : |iV-n|><J) 

= P {n+S) + 0{5e 2 ) +0(e) + (||p(0)||iProb(|iV - n| > 5)) (62) 

We claim that for 5 given by (j61|) 

Prob(|iV-n| >8) = 0{e) + <D{5e 2 ). (63) 

Indeed, since t = n(T + M) 

(n+5 \ /n-5 

J2(T + d k )<t\+ Prob I J2( T + dk)>t 
k=0 J \k=0 

= Prob ( ^=o(T + 4) - (n + 6)(T + M) ^ 5(T + M) 
\ a\/n + 5 ay/n + 8 

+ Prob ,tt(T + d k )-<^S)iT + M) > S(T + M)- 

cr\Jn — 5 (j\Jn — 5 
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We are going to show how the choice (J61|) implies 

Froh ( ZlJo(T + d k )-(n-5)(T + M) ^ 5(T + M) \ e | ^ ^ 
\ o\Jn — 5 o\Jn — 5 J 2 2 

The other half of 

Prob ( E ^° (T + dk) ~ {n + 6){T + M) C ^ T + M ) \^ | fe2 (m) 
\ aV^Td aV^+S J ~ 2 2 ' 1 ' 

is analogous. 

Depending on the size of n one has either: 

°' 8 ' <^ (67) 



a^^/n — 5 2 
or: 

"jfL. > S 4. (68) 
cr 3 \m — 5 2 

If (|67|) holds, which corresponds to large n, we use the central limit theorem with Van Beek 
rate of convergence, see [?] : 

Prob [ ZTJoiT + d k ) - (n - 6)(T + M) 6(T + M) \ < J_ [°° , 0.8p 



aVn - S a\fn - 5 / V2n J ht+m) a 3 Vn - S 

a \/n — S 



tx- 



This together with (p7jl. the inequality 



1 f°° o , e ~ a i 2 



2lX J a 



and the fact that 5 > ( T +m) V n l°g g ~ 2 implies ^ n _$ > 21ne 1 , proves (pj^j) for the case 
(|U7j). If (pjH|) holds then we apply Chebyshev inequality: 

;^(T + (i fc )-(n-(y)(T + M) 5(T + M)\ ff 2 (n - g fe^ 
aV^6 aV^6 J ~ S 2 (T + M) 2 ~ 2 ' 

where the latter inequality follows from (J68)) and 

P ) ,-6/5 



5 < 1.39 0/rT1 r — e 

From (|f)4jl. (|fi5J) and (jfifij) we get relation (jHSJ). The latter plugged into 1)62)1 proves the 
first stage. 

Finally, the three stages imply Theorem 13 . 21 provided that both e and Se 2 are dominated 
by Cmax{ne 3 ,£ 3 / 4 }, for an appropriate constant C > 0. This follows directly from e < 1 
and (pTTj) . The proof is now complete. [] 
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4 Diffusion of power in systems with discrete and con- 
tinuous spectrum 



Thusfar we have considered with systems with Hamiltonian, Hq, having only discrete spec- 
trum. We now extend our analysis to the case where Ho has both discrete and continuous 
spectrum. Continuous spectrum is associated with radiative behavior and this is manifested 
in a dissipative correction to the operator (j21j) . entering at 0(e 2 ). Therefore, the dynamics 
on time scales n ~ e~ 2 is characterized by diffusion of energy among the discrete modes 
and radiative damping due to coupling of bound modes to the "heat bath" of radiation 
modes. 

The hypotheses on the unperturbed Hamiltonian Ho are similar to those in [?]. There 
is one exception though, the singular local decay estimates are replaced by a condition 
appropriate for perturbations with continuous spectral components, see Hypothesis (H7') 
below. For convenience we list here and label all the hypotheses we use: 

(HI') Ho is self-adjoint on the Hilbert space 7i. The norm, respectively scalar product, 
on TL are denoted by || • ||, respectively (•,•). 

(H2') The spectrum of H Q is assumed to consist of an absolutely continuous part, 
cr cont (if ), with associated spectral projection, P c , spectral measure dm{^) and a discrete 
part formed by isolated eigenvalues Ai, A 2 , . . . , A m (counting multiplicity) with an orthonor- 
malized set of eigenvectors if)i, ip2, ■ ■ ■ , Vv ; i- e - for A;, j = 1, • • • , m 

Hoipk = A fc ^ fc , (ipk,ipj) = hj, 

where 5kj is the Kronecker-delta symbol. 

(H3') Local decay estimates on e~ lHot : There exist self-adjoint "weights", w_, w + , 
number r\ > 1 and a constant C such that 

(i) w + is defined on a dense subspace of TC and on which w + > cl, c > 

(ii) w_ is bounded, i.e. u>_ G C(H), such that RangeiwJ) C Domain(w + ) 

(iii) w + W- P c = P c and P c = P c W- w + on the domain of w + 
and for all / G 7i satisfying w + f G H we have 

\\w-e- iHot P c f\\ <C (t)- r ^\\w + f\\, tel. 

The hypotheses on the perturbation are similar to the ones used in the previous sections 
for discrete systems, namely: 
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(H4') (3 is a bounded self adjoint operator on 7i and satisfies \\/3\\ = 1. In addition 
we suppose that (3 is "localized", i.e. w + j3 and w + (3w + are bounded on H, respectively on 
Domain{w + ). 

(H5') go(t) G L 2 (R) is real valued, has compact support contained in [0,T] on the 
positive real line and its L 1 -norm, denoted by ||go||i is f. Therefore its Fourier transform, 
go is smooth and ||go||oo < 1- 

(H6') d ,d i, . . . are independent identically distributed random variables taking only 
nonnegative values, with finite mean, M, and such that for any / G {0, 1, . . .} and j ^ k G 
{1, 2 . . . , m} we have 

E ( e i(A ^ Afc)di ) = 

where E (•) denotes the expected value. 

Define the common characteristic (moment generating) function for the random vari- 
ables do + T,di + T, . . . 

p(0 = E ( e -^o+T)) = E ( e -*(«in-r)) = . . . (69) 

Note that p is a continuous function on M bounded by 1. Then (H6') is equivalent to 

p(X k - Xj) = 

for all j ^ k G {1,2, . . . ,m}. 

We require an additional local decay estimate: 

(H7') There exists the number r 2 > 2 such that for all / G Ti satisfying w + f G H and 
all Afc, Xj, k, j = 1, . . . ,m we have: 

|| W _e-^V(^o - A fc )so(#o - A fc )£ (A, - ff„)Pc/|| < ^p|K/||, t G R. 

Here go denotes the Fourier Transform, see JUJ), and the operators p(i?o — A)P C , #o(A — 
/Jo)Pc are defined via the spectral theorem: 

p(F -A)P c = f p(£-A)cfm(£) 

= e -Wo-\)T E ( e -i(ff -A)«) , / = 1,2,..., (70) 



£o(A-# )Pc = / 0o(A-f)dm(£) 



= / So(t)e- l(A - Ho) 'P c ^, (71) 
Jo 

where dm(£) is the absolutely continuous part of the spectral measure of Hq. 
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Remark 4.1. Conditions implying (H7') If H = — A + V(x) is a Schrddinger operator 
with potential, V(x), which decays sufficiently rapidly as x tends to infinity, then either 

E( e iA i*) = 0, I = 0,1,... and j = 1,2..., m (72) 

or 

^o(Aj-) =0, j = 1, . . . ,m (73) 

imply (H7') ; provided the mean and variance of the random variables do, d±, . . . , are finite. 
Note that {72j) is equivalent to adding the threshold, Ao = ; of the continuous spectrum to 
the set of eigenvalues {A& : k — 1,2, ... , m} for which (H6') must hold. Hypothesis \73j) 
means that the perturbation should not induce a resonant coupling between the bound states 
and the threshold generalized eigenfunction associated with Ao = 0. 

In analogy with the case of discrete spectrum, we write the solution of (J2J) in the form 

in 

(f){t,x) = ^2<lj(t)l()j(x) +P c (f)(t,x). 
3=1 

Recall that the expected power vector P(t) is defined as the column vector 
Pit) = (E(a iai (t)),E(a 2 a 2 (t)), E{a m a m {t))) . 

We denote by 

P^ = P(t), t„_!+T<t <t n 

the expected power vector after n > 1 defects (note that Pit) is constant on the above 
intervals) . 

We will show that the change in the power vector induced by each defect can be ex- 
pressed in terms of a power transmission matrix 

% = T disCtE - e 2 diag[7i, j 2 , ■ ■ ■ , 7m] 

= I — e 2 B — e 2 diag[ 7l , 72, . . . , 7m] (74) 

Recall that T disc , £ = T e = I-e 2 B, displayed in (l21H22jl . is the power transmission matrix for 
systems governed by discrete spectrum. Each damping coefficient 7 fc > 0, k = 1, 2, . . . , m 
results from the interaction between the corresponding bound state and the radiation field. 
In contrast to the results in [?], there are no contributions from bound state - bound state 
interactions mediated by the continuous spectrum; these terms cancel out by stochastic 
averaging. 
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Remark 4.2. For sufficiently small e we have: 

||T e ||i = 1 - e 2 min{7i, 72, ... , 7 m } < 1 (75) 
The damping coefficients are given by: 



Ik = lim 



9o(H - \ k )^I-\p(H -\ k -ir])\ 2 (I - p(H - A* - in))' 1 P c [(3^ k 



2 



>0, 
(76) 

for all k = 1, 2, . . . , m. Here the operators which are functions of H are defined via the 
spectral theorem and I is the identity on 7i. 

The following theorem is a generalization of our previous result on the effect of a single 
defect on the mode-power distribution, adapted to the case where the Hamiltonian has 
both discrete and continuous spectrum: 

Theorem 4.1. Consider the Schrddinger equation 

idt<f> = #o0 + g(t)/3<f>, (77) 

where g(t) is a random function, defined in terms of go(t), given by O). Assume that 
hypotheses (Hl'-H7') hold. Consider initial conditions for HJ) such that w + P c <f)o £ "H- 
Then there exists an e > such that whenever \s\ < Eo the solution of (HJ) satisfy: 

p(n+l) = % p(n) + Q ( £ 3) + Q , n = 0, 1, 2, . . . , (78) 

where the matrix T e is given in \74\) and r = min{r!,r 2 — 1} > 1. 

By applying this theorem successively we get the change over n > 1 defects: 

n-l , , » s 

(79) 



p(n) = T £ n P(0) + J2 Te (0(E 3 ) +0( 



((n-k)TY 



k=0 

Using ||T e ||i < 1 and 

00 

J>T)- < 00 

71=1 

we can conclude that the last correction term in (|79*j) is of order 0(e). 2 As for the other 
correction term we have two ways in computing its size. The first is based on ||^ fe ||i < 1, 
and gives 

-3\ 



k=0 



2 one can actually show that J2k=o 1~ E k O ^ (( n _fc)j)r J = ofnhn |e, — |^ . However, as n — > 00 the 
other correction term dominates and the result of Theorem 14.21 cannot be improved. 
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The second is based on 

n— 1 1 
fc=0 ' 

where 7 = min{7i, 72, • • • , 7m}, and gives 

n-l 
fc=0 

We have proved the following theorem: 

Theorem 4.2. Under the assumptions of Theorem \4-l\ the expected power vector after n 
defects, n = 1, 2, . . . , satisfies: 

P (n) = r £ n P(0) + (min(£7- 1 , ne 3 )) + 0(e). 
i/ere ; 7^ zs i/ie diffusion/ damping power transmission matrix given in \7$ . 
Moreover, the argument we used in the proof of Theorem 13.21 now gives 

Theorem 4.3. Under the assumptions of Theorem \4-l\ the expected power vector at a fixed 
time t, < t < 00 satisfies: 

P(t) = T £ n P(0) + 0(e 4/5 ). (80) 

Here, n is the integer part oft/(T + M), T is the common time span of the defects and M 
is the mean of the identically distributed random variables do, d±, . . . . 

The nicer form of the correction term in (fKU|) compared to (jf)9|) is due to the fact that 
min(t£ 3 , e/j) is now dominated by 0(e 4//5 ). 

In analogy with Corollary 13.11 we have, in the present context, the following limiting 
behavior: 

Corollary 4.1. Under the assumption of theorem \4-l\ the following holds: 

( P(0), ift<£" 2 

lim P(t) = { e-( B+r ) T P(0) if t = re~ 2 , (81) 

(O, if t > e~ 2 , e -> 

where B is displayed in jf^lj) and 

T = diag [71,72, ... ,7m] > 
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Proof Since T £ = I — e 2 (B + T) and B + V is self adjoint with 



J B + r>min{7 fc : k = 1, 2, . . . m} > 

we have 

{I, if n < e~ 2 

e -(B+r)r ii n = re- 2 . (82) 

ifn»£- 2 ,5^0 

This follows from writing % in the basis which diagonalizes B + T and using the fact that 
all eigenvalues of B + V are strictly positive, see the proof of Corollary 13. II 

Clearly, (JH2*|) and Theorem 14.31 imply the conclusion of the corollary. 

Note that on time scales of order 1/e 2 the dynamical system is now equivalent to: 

d T P(r) = (-B-T)P(r), 

where —B is a diffusion operator, see the discussion after relation (j51|) . while —T is a 
damping operator. 

It remains to prove Theorem 14.11 
Proof of Theorem 14.11 Consider one realization of the random variables do,di, . . . . For this 
realization the system (0) is linear, Hamiltonian and deterministic. It is well known that 
such systems have an unique solution, <ft(t), defined for all t > and continuously differen- 
tiable with respect to t. Moreover 

= WM- (83) 

We decompose the solution in its projections onto the bound states and continuous spec- 
trum of the unperturbed Hamiltonian: 

m 

<f)(t,x) = ^ajitty + P c 0(t) = (j) b {t) + d (t), (84) 
where 06 and (fid are, respectively, the bound and dispersive parts of 0: 

m 

<t>b{t) = ^%(t)V>i, 

i=i 

Mt) = Pc0(t) (85) 

and 

(<f> b (t),Mt)) = 0- (86) 
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Note that flEHJ) and flUJ) imply 

Hb{t)\\ < Hoi \\Mt)\\ < Hoi (87) 

for alH > 0. Consequently, 

k(*)| < Il0o||, (88) 

for all t > 0. 

By inserting (jHljl into (j2J) and projecting the later onto the bound states and continuous 
spectrum we get the coupled system: 

id t a k (t) = X k a k (t)+eg(t)(^(3(f) b (t))+eg(t)(^ k ,f3Mt)), (89) 
id t (j) d {t) = H 0( f> d {t) + eg(t)P c p<f> d (t) + eg(t)P c /3<f> b {t), (90) 

where k = 1, 2, . . . , m. Duhamel's principle applied to (|9U|) yields 

Mt) = e~ lHot M0) ~ ie f g(s)e- lHo ^P c (3Ms)ds - is f g{s)e- lH ^P c ^ b {s)ds. 

Jo Jo 

(91) 

In a manner analogous to the one in [?] we are going to isolate (pa in 1)91)1 . Consider the 
following two operators acting on C(R + , Domain(w + )) respectively C{R + ,H), the space 
of continuous functions on positive real numbers with values in Domain(w + ) respectively 
H: 

K + [f](t) = f g(s)w^ H ^P c (3w + f(s)ds 
Jo 

K[f](t) = f g{s)w-e- iH ^P c (3f{s)ds. 
Jo 

Then, by applying the W- operator on both sides of fpTTjl we get: 

w_Mt) = W-e-^V^O) - teK + [w„M(t) ~ ieK[<j) b \{t) . (94) 
On C(M. + ,7i) we introduce the family of norms depending on a > : 



(92) 
(93) 



and define the operator norm: 



= sup<f)°||/(f)|| (95) 
t>o 



sup \\Af\\ a . (96) 

ll/lk<i 



The local decay hypothesis (H3') together with (H4') and (H5') imply: 
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Lemma 4.1. IfO < a < T\ then there exists a constant C a such that 



\\K + \\a 5- C a 
ll-^IU — C a . 

Proof of Lemma | 4. 11 Fix a, < a < r\ and / G C(M+, Domain(w + )) such that ||/|| Q < 1. 
Then 

ft 



(*ni^ + [/](*)ii = 



[ g(s)w.e- iHo(t - s) P c (3w + f(s)ds 
Jo 



< (t) a [ \g{s)\\\w_e- iHo{t - s) P c w-\ 
Jo 



\w+Pw + \\ ■ ||/(s)||ds 



< (tyc\\w + pw + \\ 



9W 



■\\f(s)\\ds, 



where we used (H3'). Furthermore, from ||/|| a < 1 and ||w + /3iu + || bounded, we have 

-t 



(t) a \\K + [f](t)\\ < c(ty 



< c(ty 



9W 



o (t-sMsy 



■(s) a \\m\\ds 



< 



r 

c(tr Yl J 



mm(t,tj+T) 



-als 

t 

\g(s) 



{r-h<t} 



(t - s) r i(s) Q 



-als. 



By the mean value theorem 

*mm(t,i,+T) 



for some 



Hence 



We claim that 



tj < £j < min(£, tj + T). 

(tT\\K + [m)\\<c(tT Yl (t-^^y 

{r-tj<t} 



(97) 
(98) 

(99) 



(t <D a {t)- 
{r-tj<t} 

for some constant D a independent of t. This is a consequence of the fact that we are 
computing the convolution of two power-like sequences. For a more detailed proof we 
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decompose the sum into two, first running for tj < t/2 and the second for t/2 < tj < t. For 
the former we have : 

E (t-tr'Htj)- 1 * < (0 " E ^>~ Q 

{i:t,<t/2} X 7 {j:t 3 <t/2} 

^ -n 



Ms) E ^T)- a (100) 

{j:jT<*/2} 



, v — n I £ \ max(0,l— a) 



since ri > max(l,a) and tj > i,- > (j — 1)T, see (H3'), the hypotheses of this lemma, 
respectively (j97|) and (jSJ). The remaining part of the sum is treated similarly: 



{i:t/2<t J <t} {j:t/2<t J <t/2} 

't 



^ (l) E ^ T > _ri ( 101 ) 



{k:kT<t/2} 

since r\ > 1 and t — tj > kT where is such that t k +j = max{t p : t p < t}, see f!97j) and 
©• 

Now (jlOOj) and (J101J) imply (}9"§j) which replaced in (J98|) proves the required estimate for 
the operator. For the K operator the argument is completely analogous. j 

We are going to use Lemma 14. II for a = and a — r%. For Cq and C n defined in the 
Lemma, let 

C K = max {C , C ri } 

Then, for e such that Ck£ < 1, the inverse operator (/ — ieK + )~ 1 exists and it is bounded 
in the norms (j9Tlj) for a = and a = r\. Then (I94|) implies: 

W-fa(t) = (I- isK+Y 1 [w„e- iHot M0)] (t) -is (I- isK + Y l K[cf> b )(t) 

= O((t)- r ^\\w + M0)\\) -ieK[<p b } + 0{e 2 \\K\<p b }\\). (102) 

Thus we have expressed the dispersive part, <fid{t) as a functional of the bound state part, 
<p b {t). Substitution of (fTU2"J) into (JEH gives, for k = 1, 2, . . .: 

m 

fta fc (t) = -i\ k a k (t) - ieg{t) ^ (^ fc , ^Vj) 

3=1 
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- e 2 g(t)(w + ^ k ,K[cf> b }(t)) (103) 
+ £ ^(t)(O(|| W+ d (O)||(t)-'- 1 )+O( £ 2 ||K[0 6 ]||)) fc = l,2,...,m. 

In particular (|103|) implies 

n-l 

a*(t„) = e^-'^ft) + e ^ e^-^Z^do, di, ■ ■ ■ , dp), (104) 

P =i 

for all A; = 0, 1, . . . m, n > 2 and / < n. Here each constant D p depends on the realization 
of do,d\,... dp and does not depend on the realization of any other random variable. In 
addition all D p are uniformly bounded by a constant depending only on Ck above and the 
initial condition 0(0). Hence 

a k (t n ) = e-°*^-^a h (ti) + 0(e\n - l\) (105) 



for all n, I — 0, 1, 2, . . . , and k = 1,2, ... ,m. 

We multiply both sides of (|103|) with a k , then add the resulting equation to its complex 
conjugate. Then we integrate from t n to t n + T and obtain for /c = 1, 2, . . . , m 



a k a k (t n + T) - a k a k (t n ) = R 1 + R 2 + R 3 , (106) 



where 



ft n +T 

R x = -i £ Y^(ip k , (3ipj) / g(t)a k (t) aj (t)dt + c.c. (107) 

j = l J tn 

R 2 = -s 2 (w + f3ip k , g(t)a k (t)K[(f) b }(t)dt} + c.c. (108) 

R 3 = 0(e(t n )-^)+0(e 3 ). (109) 

If we neglect the R 2 and R 3 in (|10fi|) we are left with R\, which is precisely the expres- 
sion associated with the power transfer in systems with discrete spectrum; see Section El 
Moreover R 3 has norm asserted in (|78|). So, it remains to show of R 2 that 

/ rt n +T 

E (w+Pfa, / g(t)a k (t)K[(f> b ](t)dt) + c.c. 



lk PZ + 0((nT)- r ) + 0(e), (110) 



where ^ k is given by (|7fi|) and r = min{rx, r 2 — 1} > 1. 
We use integration by parts. Let 

-t 



K[(j> b ]{t)= [ g(s)e- iX ^- t - ) K[<p b ](s)ds, t n <t<t n + T. (Ill) 

Jt n +T 
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and note that K[(f> b ]{t n + T) = 0. Lemma E~T1 together with 

9(s) =g (s-t n ), t n <s<t n + T, (112) 
imply the existence of a constant C with the property: 

\\K[<j>A{t)\\< C|bo||? = C, (113) 
uniformly in t n < t < t n + T. Define 

A fc (t) = a fc (t)e lA ^"H (114) 

for k — 1, 2, . . . , m. Note that 



A k {t n ) = a k {t n ) (115) 

From ()103|) we have 

|d t A fe (s)| <C \e\ \g (s-t n )\ (116) 
for some constant C independent of s and t n < s < t n + T. Now 

t n +T ftn+T 

g(t)a k (t)K[<f) b }(t)dt = / A k (t)d t K[(j) b \(t)dt 

rt n +T 

= -a k {t n )K[(j) b \(t n ) - / $A fc (*)#[&](*)d* (117) 



a*(0 / ^ + %(t)e a *<'-^lir[&](t)dt + 0(e) 



To further rewrite ()117j) we note that for t n < t < t n + T 



K[<Pb](t) = Y. a 3 (s)g(s)w.e- tHo ^P c Pij j ds 

j=i n 

n-l m ^t,+T 

+ a j {s)g{s)w^ lHo(t - s ^ c f3^ j ds (118) 



l=0 j = l Jt l 

An integration by parts similar to the one above and use of ()112|) leads to: 



a j {s)g{s)w^e- lHoi - t - s) Y t c ^ j ds = 

a]itl) /" +T ^ s ) e ~ iAj(S ~" )w - e " // ° ( ^ S)p ^^ + ° ( {t _ t "_ T y, 
a^W-Uh ~ H )e-^ t ~^P c (3^ + O ( (t _^_ T)n ) , (119) 
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and 

h 

a j (s)g(s)w^e- iHo ^P c ^i) j ds 



= aj(t n ) [ g(s)e" iX ^ s - t "~>w_e- iHo{t - s) P c (3^ j ds + 0{e). (120) 
Jt„ 

By plugging (jll9H120|) in ()118|) we get 

m „ t 

3=1 Jtn 

m n— 1 

+ E E a A t i) w -9o{\ - ffo)e- iflo(i - i!) P c ^ + O(e), (121) 

j=l (=0 

where to estimate the error we used the fact that the series ~~ h ~ T)~ ri is convergent 
and uniformly bounded in t. 

We now substitute (J121)) into the right hand side of (|117j) and obtain 

g{t)a k {t)K[(j> h ]{t)dt = 0{e) 

m r t n +T 

+ ^cikit^a^tn) / gity^-^ / g(s)e- tX ^ s - tn) w^e- tHo{t - s) P c f3^ j dsdt 

j — l ** tfi J t n 

m n— 1 

+ E ^Mtn)"j(ti)v-9o( H o ~ h)g (Xj - Fo)e~ !ff,,(i "" i,) P c W r (122) 

3=1 1=0 

Based on (jl(J4j) we can replace a k (t n )aj(ti) in ()122|) with 

o fc (t n )aj(*0 = e iAfc ^-*^a fc (^)a i (^) + error(Z,j) (123) 

n-l 

errar(l,j) = eJ2^ k{tn ' tp) D P W-g Q (H - X k )g (X, - tfoK^^Pc/^- 
P =i 

Taking into account that t n — t n -\ = d n + T and the fact that £ n _i — ti, D p , I < p < n — 1 
do not depend on d n , the expected value of the error can be rewritten as 

K(error(l,j)) = 
n-l 

= eJ2 E (^-e i{Afc " Ho)(t ' l "*"- l) e afe(i "- 1 -* p) J D^o(^o - A fe )^ (A J - #o)e~ iHo(t "- 1_ ' !) Pc/^) 
P =i 

n-l 

= eJ2w-p(H - X k )E (e iX ^-^D p g Q (H - X k )g (X 3 - H )e~ mo(t ^-^P c f3^) 
P =i 

n-l 

= eJ2 E (e iX ^-^D p w-p{H Q - X k )g (H - X k )g (Xj - H ) e - iHo ^-^P c f3^) (124) 



P =i 
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By applying the Ji norm to (J 124)) . commuting the norm with both summation and expected 
value and using (H7') we get: 

\\E(error(l j3 ))\\ < \e\ ^ ~ ^ < C \e\((n - l)T)^. (125) 

Since r2 > 2 the summation over I and j of all the errors will have an 0(e) size. By this 
argument ()122j) becomes: 

-t n +T 

g(t)a k (t)K[<p b ](t)dt} + c.c.) = 



E (Jw + M h ,£ g(t)a k (t)K[(f> b }(t)dt^ + c.c)j = ^E(a fc (* n )o i (t„)) 

■ (w+Pipk, g{t)e lXk{t " tn) j g(s)e- iX ^ s - tn) w.e- iHo{t ~ s) P c ^ j dsdt) + c.c. 

m n— 1 

j=i ;=o 

■ E ((w + M k , w.g (H - A fc )po(Ai - i^e^""^-''^^^)) + c.c. 

+ 0(e). (126) 

But (H6') and the technique used to prove (J4*2"j) imply 

"<*(«)-,(«)) = { f 

Moreover, an argument similar to the one we used in (jl23H125j) allows us to replace by 
in flZBj) and incur an 0(e) total error. Then, (|12fij) becomes 

E (Jw + Pi> k , £ g(t)a h (t)K[<j> b ](t)dt\ + c.c^j = 
= P^ (w + (3i/j k} w_g (H - A fc )^o(A fc - H )P c (3ip k } 

+ Pi n) (w+Mh, w-g (H - \ k )g (\ k - H )E fa e i(A k -^o)(*„-t I )p c j + c . c . 

+ O(e). (127) 
We claim that 

7^ d = (w + ptp k , w^g (H - \ k )g (\ k - H )P c p?p k ) 

+ (w + M k , w„g (H - \ k )g (\k - H )E fa e ^ k -H )(t n - tl ) p \ ^\ + c c 

= lk + 0{(nT) 1 -^) (128) 
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where j k is given in (J76|) . (|128|) replaced in (J127)) gives (jllO)) which finishes the proof of 
this Theorem. 

To prove (jl28j) we first find a simpler expression for the expected value operator involved. 
Since {dj}j>o are independent, identically distributed with common characteristic function, 
p(£), using the definition of t n , n > 0, see (jSJ) and the spectral resolution of the operator 
H , see flTUj) . we have: 



J Ccont (Ho) 



(Ho) 



& cont (Ho) 

n-l 



n IL J. 

= / IjE(e i < A *-flte +r >) dm(£) 

= / p n ~'(£ - A fc ) dm(0 = p n ~\H - X k )P c . (129) 

J 0"cont (-Wo) 

Hence 

w_£o(#o - A fe )^ (A fc - F )E ^ e ^-*°)(*«-*0p c j (3 

n 

= w_g (H Q - X k )g (X k - H Q )^(P{H Q - X k )P c f3. (130) 

But each operator term in (|13Uj) has its 7^— norm dominated by: 

\\w_g (H - X k )g (X k - H )p\H - X k )P c (3\\ = 
= \\w-p(H - X k )g (H - X k )g (X k - H )E(e-^ H °- x ^-^P c (3)\\ 

^ ,+ c , v , \M\\<(u-m- r2 - 

Now r 2 > 2 implies that the sequence 1/ (jT) 7 " 2 is summable, and, by the dominant conver- 
gence theorem, there exists: 

7 fc = (w + (3i[) k , W-g (H - X k )g (X k - H )P c P^k) 

00 

+ ( w +^k, w^g (H - X k )g (X k - H )fP{HQ - X k )P c f3ip k ) + c.c. 

3=1 
= it- 
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Moreover 

oo 

\lk~lk\ = ^2 ( W +P^ k ' w -9o(H - \ k )g (\ k - H )p J (H - \ k P c f3ip k ) + c.c. 

j=n+l 



(131) 



]=n 



Consider now, for rj > 0, 

H = (w + f3?p k , w_g (H - \ k )g (\ k - H )P c pip k ) 

oo 

+ ( w +^k, w^g (H - \ k )g (\ k - H )p j (H - A fe - ir])P c f3ip k ) + c.c. (132) 



On one hand 

p>(H - A fc - ir))P c = E(e- r ' {t >- to) e- i{Ha - Xk)(t >- to) P c ) 
and, by the dominant convergence theorem, for all j • > 1 

lW(tf - X k - ^)P C = p>'(#o - A fc )P c . 

r?\0 



(133) 



On the other hand the series ()132|) is dominated uniformly in rj by a summable series, 
because: 

\\w-g (H - h)go(h - Hq)p>(H q - A fc - ^)P C /3|| 

/ / dudsg (s + u)g (u)E ( e -'^-*° Ve- l(H °- Afc)( ^-*°- s) P c /3) 
Jo Jo 



< 



Ce~™ T 



{t. -t -T) 



-IbollilK^II < (O'-iW 



■ri 



Here we used (H3'), \\go\\i = 1 and bounded. Therefore, by the Weierstrass crite- 

rion: 

lim7 n k =% (134) 

7?\0 

In addition ()133j) implies 

\\p(H -X k -iri)P c \\ < E(e- ,?( * 1 -* o) ||e- i(Ho - Afc)( * 1 -* o) P c ||) 

< e-" T <l. 

This makes (I — p(Hq — X k — ii]))P c invertible and given by the Neumann series: 

oo 

(I - p(H - Afe - ^))- 1 P c = £y (ff - Afe - i7j)P c . (135) 

j=0 
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Plugging (|135|) in (|132|) we have 



Ik 



((3ip k ,g (H Q - \k)go(h - H )P c Pi) k ) 



+ (Pi>k, 9o(Ho - Afc)<?o(Afc - H q )p(Hq - \ k - irj)(J - p(H - X k - irj)) ^c/fy^) + c.c. 
A simple inner product manipulation shows that: 



72 



9o 



Hence 



{H - X k )y/l-\p(H -\ k -irj)\ 2 (I - p(H - A fc - ir,))" 1 P c [/3^ 



7 fc = lirn^ = 7fe , 

r/\0 



(136) 



see also (fH^jl and ifTfij) . 

Finally, (|136j) and p31j) give the claim (|128|l . The theorem is now completely proven. 



5 Comparison to stochastic approach 

In this section we want to compare our results with the stochastic approach in [?, ?, ?, ?]. 
We view the results of this paper and those discussed in this section as complementary. 
The results of this paper apply to the situation when a known localized "defect", g , is 
randomly distributed in a manner which achieves averaged diffusive effect. The results of 
Papanicolaou et. al. apply to a random medium, which is unknown and with assumptions 
about their distribution. One of the key technical assumptions in this latter work is that 
the expected value of the randomness, at any time, is zero, i.e. in our notation E(g(t)) = 0. 
In the results of this paper, we allow for E,(g(t)) to vary with t. Indeed, for our train of 
pulses (see (jlj) and figure |TJ) E(g(t)) = and implies g (t) = 0, so unless we have the 
go = 0, K(g(t)) is generally different from zero and time-dependent. On the other hand, 
our hypothesis (H4) has no corresponding restriction in Papanicolaou et. a/.'s theory. 

Another important difference is that our result applies on time scales even larger than 
1/e 2 , where e is the size of the randomness while the other results apply only on time scales 
up to 1/e 2 . However, it appears that there is a striking similarity between the two results 
on 1/e 2 time scales. The train of pulses we analyzed is closest to the stochastic process 
described in [?, Section 2] where both its values in the epochs [0, d + T], [d + T,d + T + 
d% + T] . . . , and the epochs are now dependent on the realizations of the same random 

variables, do, d\, However, if we assume that the radiation modes are not present, the 

dynamical system we investigate, ()89II90|) . is the one in [?, Section 4], see also [?]. These 
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prevents us to use the formulas in the above papers. Nevertheless, we are going construct 
another stochastic perturbation, in the spirit of the one in [?, Section 4], for which we can 
compute the expected power evolution using both theories. We find that the two results 
coincide but keep in mind that while the example satisfies all our hypothesis it does not 
satisfy one of theirs, see below. 

In addition to (H1)-(H4), suppose the random variables do, d\, di, • • • , can only take 
values in the interval [0, d] (this is clearly satisfied by the random variables constructed in 
the previous section for finitely many modes) and denote by //(£) the measure induced by 
their distribution. Consider the positive real axis partitioned into "epochs": [0, T+d],[T+ 
d, 2(T + d)], . . . of length T + d. In each epoch a defect is placed at a distance dj from the 
starting point of the j th epoch. Specifically, the first defect is placed at a distance do from 
t = 0, the second at a distance d\ from t = T + d, .... Here do, d±, d 2 , ■ ■ ■ are realizations 
of the random variables do, d\, d 2 , • • •■ That is, we will now consider equation (|12|) with the 
perturbation given by: 



g(t) = g {t - d ) + g (t — (T + d) — d{) + g {t - 2(T + d) - d 2 ) + . . . . 



(137) 



see figure |21 



9(<) 



90 (t - d ) 



90(< - (T + d) - di) 



90 (i - 2(T + d) - da) 




T + d 



di 




2(T + d) 




d-2 



3(T + d) 



Figure 2: Another train of short lived perturbations 

Our result, Theorem 13.11 applies without any modifications since before each perturba- 
tion we have: 

E (a k a,j(l(T + d) + dt)) = E (a k aj(l(T + d))) E (e^*^ 1 ) = 0, 

if k ^ j. As for Theorem 13.21 its proof is much simplified and the error estimate improved 
because we now know how many complete defects are going to appear up until the chosen 
time "i" , namely n = \ t/(T + d)\ . The expected power at time t can differ from the one at 
time n(T + d) by no more than the size of the perturbation, e, since after each experiment 
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only a part or a single full defect can occur in between this time slots. Hence: 

P(t) = P { n) + 0(e) = TJ\P(0) + 0{ne 3 ) + 0(e), (138) 

where the integer n is such that n(T+d) <t< (n+l)(T+d). To get closer to Papanicolaou 
et. a/.'s results, suppose 

t = t/e 2 , t > is fixed 
and pass to the limit e \ 0. We get 



limP(t) = limT £ L(J+d,e J P(0) = lim(I - e 2 B) Li^J P(0) = e TB P(0) (139) 



)i(T+d) 



e\0 e\0 

where B is given in (1221) and 



e\0 



B = -W^) B ' (140) 

Let us now apply Papanicolaou et. al. result to the above example. Note that the 
manner in which the perturbation is constructed makes the example very close to that in 
[?, Section 4]. But since the stochastic process is not piecewise constant, one has to rely 
on more general form of their results such as [?, Remark 2 in Section 2]. The ODE system 
for the amplitude vector, a(t) = (ai(t), a 2 (t), . . .)', is: 



d t a(t) = Aa(t) — ieg(t)aa(t), 
a(0) = a(0), 



(141) 



where 



A = -idiag[Ai, A 2 , ■ • .], 

a = ((ipk,^))^^] 

see also (JT7J). This is a special case of system (2.27) in [?] with 

M = -ig(t)a. 
Note that hypothesis (2.28) in [?] translates into 

= E(g(t))= I g (tf - s)dp(s), for all i > 
Jo 

where 

t' = t - (T + d) 



(142) 



t 



T + d 
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which generally implies the trivial case g = 0. Hence for nontrivial examples Papanicolaou 
et. al.'s theory is not rigorously applicable. We are going to replace (1142)1 with a milder 
one: 



1 



lim - / E(g(t))dt = 0, (143) 

t/oo t J 

i.e. the time average of the expected value of the perturbation is zero, under which we 
can formally derive closed coupled power equations. While the results of the stochastic 
approach do not apply to this example because E(g(t)) ^ 0, our results do apply and it is 
reasonable to conjecture that there is an extension of the stochastic approach to this case. 
In the special case, where g is given in (j!37|) . the condition (J143)) reduces to <?o(0) = 0, 

Let us compute their equation for the evolution of the powers, i.e. we prove that their 
system of equations for the product of amplitudes: 

d T E(a <g> a) = VE(a®a), 

where 

t = — is fixed as e \ 0. 
gives a closed equation in powers, i.e. 

V vqm , =0, if q ± q' (144) 

and consequently for the powers P(r) = diagE(a ® a(r)) we have 

d T P(r) = VP{t), (145) 

with 

Vpg = V mm . (146) 

The main point is that V coincides with B in our result (J139)) : see also ()140)) and ()22)). Thus 
the two results agree on time scales of order 1/e 2 . 

For the formula of V we only have to replace M in [?, equation (2.35)] by its complex 
conjugate M whenever it applies on the right part of the tensor product, i.e. 



1 ft+to PS 

V = lim - / / E (M(s)M(o-) ® / + M{s) ® M(o)) dads 

+ lim i [ I E (M(a) ® M{s) + 1® M{s)M{a)) dads, (147) 



to Jto 



M pq = -ie iA ^g(t)a pq ; 
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see also [?, relation (2.32)]. It will be clear from the argument below that the limit in (|147|) 
does not depend on t (note that this in in fact a requirement for the validity of the theory) 
so we are going to work with t = 0. Although the computation of V has been done in [?] 
(denoted there by V) and then summarized in [?, Section 3] we are not able to use them 
because they relied on the stationarity of the process, see [?, relation (2.2)] which is not 
satisfied by our example. Nevertheless we have component wise: 



V pqm , = -^Va^, lim ~ / f e iAprS e iAr « a E(g(s)g(a))dads 
r t y°° 1 Jo Jo 

+ a pq a q , p lim - f [ e iAp " s e iA ^'^E(g(s)g(a))dads 
t/co t J J 

+ a pq >a qp lim - [ [ e lA ™ s e iA ™' a E(g(s)g(a))dads (148) 

t/oo t Jq Jq 

- 5 pq a rp a q > r lim - / / e iArpS e iA «' ra E(g(s)g(a))dcrds, 

t/oo t Jq Jq 

where we have used ctkj = &jk due to the self-adjointness of (3 in a^j = (ipk, fiipj) and the 
fact that g(t) is real valued. Thus it is sufficient to compute 



t ps 

iAkjsAAAicr 



JO 



e tA ^ s e tA ^E(g(s)g(cr))dads. 



Let us fix t and suppose for the moment that t = n(T + d). Then 



t rs 

iA ki s iAna 



JO 



e l ^ s e l ^ lf7 E(g(s)g(a))dads 

n-l n(m+T)(T+d) ps 



"■~ j - i>(m+L){i +a) ps 

V / / e iA ^ s e iA ^E(g(s)g(a))dads 

m=Q Jm(T+d) Jo 

= / e iA ^ s e iA ^E(g(s)g(a))dads (149) 

m=0 J™{T+d) Jm(T+d) 

p(m+l)(T+d) rm(T+d) 

+ V / e iAk > s E(g(s)) / e iA ^E(g(a))dads, 

n Jm(T+d) Jo 



where we have used the fact that the random variable g(s) and g(a) are independent unless 
s and a are in between the same epochs. Now 

/ e lA ^E(g(a))da = Y) / e lA ^E(g (a - r(T + d) - d r+1 ))da 

Jo r=Q Jr{T+d) 



41 



m -! rT+d pa 

V e l ^ ir{T+d) / e iA »° / p (o- - s)dfi(s)da 

r=0 J° ^° 

m-1 /.d pT+d-s 

y^ e ^T + d) / e ,A, iS / e lA ^^ ( ( x)rf ( T^( S ) 



r=0 
m—l 



lit— i „ d „j 

yyA, iK T +( *) / e ^ rf//(s) / e ^^ o((T)AT 

m—l 

^- e iA J ,r-(T+d) E ^iA J ,d r+1 ^^_ A ^) = Q) 



r=0 

where we used supp g C [0,T], E(e iA ^ ldr+1 ) = if j ^ I, see (H4), and the fact that 
9o(0) = 0. 

The only nonzero terms left in (J149)) are of the form: 

p(m+l)(T+d) ps 

/ / e iA ^ s e iA ^E(g(s)g(a))dads 

Jm{T+d) Jm(T+d) 

rT+d rs 



pl+a ps 

= ^ kl m { T+d) / / e ^ kj s e i^^ g ^ s _ dm)go{(J _ dm))dads 

JO JO 
pT+d pT+d pd 
= e ^kMT+d) / e iA fci s e iA jiCT / go ( s - - Qdv{t)dsd<T 

Jo J a Jo 

= el A kim ( T+ d) f\iA^ [ -f+^t e **° go ( <T ) f +d ? e lA *> s g (s)dsdad^). 

Jo J J a 

Now, the upper limit of the integrals with respect to s and a can be replaced by oo without 
modifying their values since supp go C [0, T] and £ G [0,d]. Hence they do not depend on 
£ and by computing the integral with respect to the measure d/i(£) first we get: 

nd 

/ e iAfci5 d/i(0 = E(e iAkldm ) = 5 kl 
Jo 

Knowing that k — I in order to get a non zero result, we can now compute the integrals 
with respect to s and a using (|37jl without the complex conjugate part. In conclusion we 
have 

J J e lA ^ s e iA ^E(g{s)g{a))dads = 

' [_^j^( w ^ + ip.v.£J^,) +0(1) 

where the correction is needed for t ^ n(T + d). Consequently 
lim / / e iA ^ s e lA ^E{g(s)g{a))dcxds 



t/oo 



JO 
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- ^{k-^™C$$»Y (l50) 

Replacing f)150|) in the formula (jl48j) it is easy to see that V pqtPq > = unless q — q' and 
in the later case the first and the fourth terms in (jl48j) are complex conjugate which is true 
for the second and third terms also. Simple arithmetic leads to 

where B is given by (jl40}) and ()22j) . 

In conclusion, on time scales of order 1/e 2 our results for the example in this section 
coincides with the one obtained by Papanicolaou et. al. in the series of papers [?, ?, ?, ?]. 
As mentioned earlier, although our result applies directly, the stochastic approach requires 
the E(g{t)) = for all t. 

6 Appendix: Properties of the power transmission 
matrix 

In this section we prove the properties of the matrix (linear operator) B we used in Corol- 
laries 13.11 and 13.21 Recall that B is given by (|2~2jl and is irreducible, see the discussion 
before Corollary 13.21 We note that (J22|) implies in particular that: 

1. all the components of B are real; 

2. bu > for all % = 1, 2, ... ; 

3. bij < for all i, j, i ^ j; 

4 - J2j bij = or equivalently bu = - Y^jj^i bij for all i = 1, 2, 

Lemma 6.1. If the dimension of B is finite, say m, then B is a nonnegative, self adjoint 
matrix having as a simple eigenvalue with corresponding normalized eigenvector: 

r = -L(l,l,...,l)' 
v m 



Proof. Since all components of B are real, self adjointness is equivalent to 

bjk = b kj , Vj, k, j ^ k. 

From (|27|) we have 

I I 2 — l — I 2 — I I 2 

l^j'fcl l^fcjl l^fcj'l 5 
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where a denotes the complex conjugate of the complex number a. 

Now because % is the Fourier transform of a real valued function, see and (H3), 
and because A k j = —Aj k , see (pfi|) . we have 



l^o(-Aj fc )| = \go(A kj )\ = g (-A kj ) = |^ (-A fei )| . 

Hence, for all j ^ k 

b jk = - \a jk \ 2 |^ (-A jfe )| 2 = \a k j\ 2 \g (-A kj ) | 2 = b kj , (151) 

rendering B self adjoint 3 . 

In order to prove that B is nonnegative, consider an arbitrary vector X = (Xi, X 2 , . . . , X m )' 
and let X* = (X\, X2, • • • , X m ) denote its adjoint. Then 

m m 

X*BX = bijXiXj = bn\Xi\ 2 + bijXiXj 

= - £ b ij \X i \ 2 + b nXiXj (152) 
= K\ • \ x i~ x i\ 2 ^ °» 

where we used properties 3. and 4. above. Hence i? is nonnegative. 
Now, if Y = Br then 

by property 4. Consequently is an eigenvalue for S with corresponding eigenvector r . 

To prove that is a simple eigenvalue we use the irreducibility of B. On the set of 
components {1,2,..., m} of vectors in C m we define the following relation: 

Definition 6.1. We say that components i and i are always coupled to zeroth order. 
We say that components i,j are coupled to first order if bij 7^ 0. 

We say that components i,j are coupled to n th , n > 2 order if there exists a sequence 
of components kx, k 2 , ■ ■ ■ , fc n _i, such that the pairs 1, k\\ fci, k 2 ; ■ ■ ■ ; k n -i,j] are all coupled 
to first order. 

We say that components i,j are coupled if they are coupled to any order. 



3 Identity i|151|) does not rely upon B having a finite dimension. Therefore it is valid even when B has 
infinite dimension. 



44 



It is easy to show that "to be coupled" is an equivalence relation on the set of compo- 
nents {1,2,..., m}. Hence it induces a partition of the components. 

Claim 1. If B is irreducible the above partition is trivial. 

Indeed, if we assume contrary the partition is formed by at least two proper subsets of 
the set of components {1,2,..., m}. By a reordering of the components, i.e. a reordering 
of the standard basis vectors in C m , we can assume assume that the partition is formed by: 

{1,2,..., mi}, {mi + 1, mi + 2, . . . , m 2 }, . . . 

Then 6^ = whenever i,j fall in different subsets of the partition, otherwise they would 
be coupled. Consequently, B has the form: 

£? = diag [fii,^,...], 

where B\ is a m\ x mi matrix, B 2 is a m 2 x m 2 matrix, etc. But these contradict the 
irreducibility of B, see also the discussion before Corollary 13.11 

Claim 2. If X = (Xi, X 2 , . . . , X m )' is a zero eigenvector for B and i,j are coupled then 
Xi = Xj. 

Indeed, X*BX = because BX = 0, and ()152j) implies 

\*Hj\'\Xi-Xj\ 2 = 0- (153) 

i,j,i<j 

If i,j are coupled to the first order then by definition 6^ ^ and we must have Xi = Xj 
in order for ()153j) to hold. By induction on the order of coupling one obtains the result of 
the claim. 

Finally, Claim 1 and the irreducibility of B imply that all components are coupled. 
Then Claim 2 implies that all components of a zero eigenvector must be equal. Hence all 
zero eigenvectors are parallel to r . Since B is self adjoint this means that is a simple 
eigenvalue. \ 

Lemma 6.2. If B is infinite dimensional, then B is a bounded linear operator on i 1 with 
— 2. In addition, for \e\ < 1, the operator T e = I — e 2 B transforms positive vectors 
(i.e. vectors with all components positive) into positive vectors and conserves their i 1 norm. 
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Proof. It is well known that B = (bij)\<i t k<oo is a bounded linear operator on i 1 iff there 
exists a constant C > such that: 



J>d <C, Vj = l,2, 



(154) 



i=l 



In this case || i? ||i < C. We are going to show that for B given by (|22|) we can choose C = 2 
in (fT541 . 

Indeed, let us fix an arbitrary j G {1,2,...} and consider the j th vector in the standard 
basis of i 1 : 



X — (Xi, x 2 , 



iii^j 

1 if % = j 



Let 



A 



[a 



ij ) l<i j'<oo 



aij = (Xij g (-Aij) — (ipi, (3ipj) / g (t)e i{x '- x ^dt. 
By a direct calculation we have 



EM = E 



i=i 



i=l 

oo 



fc=l 



Xi dikdkpXp — dikXk o-ikXk 

k,p=l k=l 

= Y,Wi(A-AX) i -{AX).{AX) i \ 

i=l 

oo oo 
1=1 8=1 

Clearly lef 2 , ||X|| 2 = 1. We are going to prove below that: 
Claim 3. A is a bounded operator on £ 2 with ||A|| 2 < 1. 
Hence 



£|(ax)< 



2 = ii^m < ii^iii = i, 



i=i 



while using Cauchy-Buniakowski-Schwarz inequality we have: 

oo / oo \ 1 / 2 / oo 

£1^1- \{A.AX)t\ < £iXf [Y.U-AX), 



1/2 



i=l 



,j=l / \i=l 

\X\\ 2 ■ \\A-AX\\ 2 < \\A\\l\\X\ 



< 1. 



(155) 



(156) 



(157) 



(15* 



(159) 
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By plugging in (|158|) and (|159|) into (|157|) we get 

oo 

EN<2 



(160) 



i=i 



and, since j was arbitrary, (j 154)1 holds with C = 2. Consequently, B is a bounded linear 
operator on C 1 with norm ||-B||i < 2. 
Consider now 

T £ = I- e 2 B, T £ = {ti 

Then for i ^ j, 

Note that by ()151j) we also have: 



l 3 ' l<i,j<oo 



ta = -e 2 h A > 0. 



tji tiji Vz, J. 



On the other hand 

t H = l- e 2 bu = 1 - e 2 E IM> 
where we used properties 3. and 4. above. Moreover 



EN = K\ + E \ b a\ 

3=1 jjj^i 



■ E hi i 



+ Eni = 2 E>;I- 



Using now ()16(J|) and (|151|) we have 

oo 1 oo 

En = ^En = ^En<i- 

Hence 



jj^i 3=1 



3=1 



We also have: 



t H = 1 - e 2 E IM > 1 - ^ 2 > 0, if |e| < 1. 



E = *« + E = 1 - £2 E 1^1 + e2 E i^'i = x > 



and by <fT77T|> 



oo oo 



(161) 



(162) 



(163) 



i=i i=i 



47 



Now let 



x = (x 1: x 2 , ...ye e\ x j > o vj = 1,2, . . . 



Then 



CW; = E^ X i>0 
since all terms in the sum are nonnegative with at least one being strictly positive. Moreover 



oo oo 



oo oo 



\TsX\u = =EE^ = E^E^ = E^- = W 1 ' 

i=i j=i j=i i=i j=i 



i=i 



where we exchanged the order of summation because we are dealing with convergent series 
with nonnegative terms and we also used (| 163ft . 

The Lemma is now finished provided we prove Claim 3. Let 



X — (X 1 , x 2 , 
be arbitrary and denote by 



\X\\n = 1 



Y(t) = ^e^Xtfj, Y(t) E H, \\Y(t)\\ = 1. 



Then 



\X*AX\ 



E a jk X j X k 

j,k=l 



9o(t)(ET=i elXjtX ^vEtie iXkt X k PiJ k )dt 



< 



oo 
oo 



< 



g (t)\-\{Y(t),f3Y(t))\dt 

\g {t)\dt= \\/3\\ n \\g \\ L i = 1 



H 



g (t)(Y(t),(3Y(t)}dt 

\g (t)\.\\P\\ n \\Y(t)\\ 2 dt 
where, at the very end, we used (H2) and (H3). j 

Lemma 6.3. If B is infinite dimensional, then B is a bounded, linear, self adjoint, nonneg- 
ative operator on i 2 with spectral radius less or equal to 2. Moreover, is not an eigenvalue 
for B. 

Proof. Because of ()151j) B is symmetric on £ 2 . Consider the 2-form induced by B on I 2 : 



oo oo 



X*BX = bjjXjXj = bu\Xi\ 2 + bjjXjXj 

i,j=l 4=1 i=l 



E m^i 2 + E E%- v <- v * = E E M(W - x i x i)> ( 164 ) 
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where we used properties 3. and 4. above. Hence 

oooo / I v 12 I v 12 \ oo oo 

X * BX = EEn(W+ ' 1 2 m )= 2 EEn-w 2 

oo 

< 2sup(E^|%|)El^| 2 ^ 2 H X ll2' 

where we used (jl62j) . So the 2- form induced by B is bounded. Since B is symmetric this 
implies that B is a self adjoint, bounded operator on £ 2 with ||-B||2 < 2. Therefore its 
spectral radius is less or equal to 2. 
Moreover, (jl64J) implies 

X*BX = \K\ ■ \ X i ~ X i\ 2 ^ °- 

i<j 

On one hand this shows that B is nonnegative. On the other hand, together with obvious 
generalizations of Claims 1 and 2 in Lemma 16.11 for the case of infinitely but countable 
components, it shows that if a zero eigenvector for the irreducible operator B exists then 
the eigenvector should have all components equal. However such a vector is not in £ 2 unless 
it is trivial. Therefore is not an eigenvalue for B. [ 
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